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A new approach to phenomenological description of the charged particle multiplicity distributions 
in pp/pp and e + e~ collisions is presented. The observed features of the data are interpreted on the 
basis of stochastic-physical ideas of multiple production. Besides the processes of parton immigration 
and absorption, two and three patron incremental and decremental recombinations are considered. 
The complex behaviour of the multiplicity distributions at different energies is described by four 
parametric generalized hypergeometric distribution (GHD). Application of the proposed GHD to 
data measured by the CMS, ALICE, and ATLAS Collaborations suggests that soft multiparton 
recombination processes can manifest itself significantly in the structure of multiplicity distribution 
in pp interactions at very high energies. 
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I. INTRODUCTION 



The multiple production of particles in high energy collisions received a lot of attention during many years. A 
renewal interest in this topic originated recently with the beginning of operation of the Large Hadron Collider (LHC) 
at CERN. Among first results obtained at the LHC were multiplicity measurements in proton-proton collisions. The 
CMS, ALICE, and ATLAS Collaborations provided valuable and precise data on multiplicity distributions (MD) 
of the charged hadrons in the new super high energy domain. The study of multiparticle production can give 
important information about parton processes and transitions of the colored partons to the colorless hadrons. The 
complex phenomena that influence MD are however very hard to describe in detail, because we cannot formulate the 
multiparton process in terms of its QCD theory. Most of the particles which contribute to the multiplicity buildup 
are soft particles and a perturbative approximation is inapplicable here. Second, hadronization should be somehow 
taken into account which is far of being understood properly. The natural and probably most economical approach 
is to look for empirical relations and regularities. 

The study of particle production as a function of multiplicity has revealed a very popular Koba-Nielsen-Olesen 
(KNO) scaling [1]. In pp/pp collisions, the scaling in the full phase-space holds up to the highest energy of the 
CERN Intersecting Storage Rings (ISR) but it is clearly violated Q in the energy region of the CERN Super Proton 
Synchrotron (SPS) collider and beyond. The strong violation of KNO scaling was observed also at the energy 
y/s = 7 TeV in the limited pseudorapidity intervals (\ri\ < 2.4) though for small pseudorapidity windows (|fy| < 0.5) 
the scaling is approximately valid [3|. A related phenomenon is the so-called negative binomial regularity, which is 
the occurrence of the negative binomial distribution (NBD) in different interactions over a wide range of the collision 
energies. The UA5 Collaboration showed that MD in the non-single-diffractive (NSD) pp/pp collisions can be described 
by NBD up to the energy ^/s = 546 GeV both in the full phase-space [J] and in symmetric pseudorapidity windows 
Q . Analysis of data on MD measured by the ALICE Collaboration indicates 0, 0| that NBD describes the data in 
the small pseudorapidity window \rj\ < 0.5 up to the energy y/s = 2360 GeV. Besides pp/pp collisions, NBD has been 
applied to various systems including e + e~ annihilations. Much effort has been made to explain the negative binomial 
form of MD observed in many situations, however, its physical origin has not been fully understood [8j. 

The shape of MD of particles produced in hadron-hadron collisions at high energies is quite different. The full 
phase-space data on charged particle multiplicities obtained from the NSD events in pp collisions at y/s = 900 GeV 
Q indicated that besides KNO, the negative binomial regularity is violated as well. The measurements of MD at the 
energy — 1800 GeV by the E735 Collaboration [13] at the Fermilab Tevatron showed even stronger deviation of 
data from NBD. Though with large mutual discrepancies at high multiplicities, both data demonstrate a narrow peak 
at the maximum and some structure around n ~ 2 < n >. Despite its correct qualitative behaviour, the single NBD 
is not sufficient to describe the experimental data. A systematic study of the complex form of MD was performed 
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in the framework of two component model [ll| and its three component modification [12j. There exist also other 
approaches to explanation of the observed features of particle production at high energies (for a review see Ref. [H). 
New interest in this field is motivated by the recent results of the multiplicity measurements in pp collisions at the 
LHC. The data obtained by the CMS 0, ALICE pj], and ATLAS [H| Collaborations show similar structure in MD 
of charged particles produced in the limited windows in pseudorapidity. The measurements allow to study evolution 
of the distinct peak at maximum and the broad shoulder at large n both with the collision energy and pseudorapidity. 

In this paper we propose an alternative phenomenological approach to the description of MD in pp/pp and e + e~ 
collisions. Using the same concept for hadron and lepton interactions, we try to account for the structure in data which 
emerges in the pp/pp interactions at high energies. The proposed four parametric representation of MD is motivated 
by a scenario of the parton cascading. The considerations behind are based on a premise that dynamics of the particle 
production as manifested on the level of multiplicities can be described in terms of a stochastic cascading with specific 
types of the underlying processes. The construction includes the processes of parton immigration and absorption 
together with two and three parton incremental and decremental recombination. The recombination processes are 
assumed in the final stage of the parton evolution during the color neutralization. The transition to the colorless 
hadrons is considered in a stationary regime at breakdown of the confinement and onset of hadronization. 



II. RECURRENCE RELATION BETWEEN MULTIPLICITIES N AND N+l 

The number of particles created in high energy collisions varies from event to event. The distributions of probabilities 
P n of occurrence of the multiplicity n provide sensitive means to probe dynamics of the interaction. Besides general 
characteristics of particle production the distributions contain information about multiparticle correlations in an 
integrated form. The correlations of the final particles reflect features of hadronization mechanism and properties 
of the QCD parton evolution just before hadronization. The data on MD of charged particles from the high energy 
pp/pp and e + e+ interactions allow us to study the processes underlying the multiple production and its correlation 
structure within a phenomenological framework of parton cascading. The extraction of information on these processes 
requires the elementary charge conservation to be taken into account. In our approach we investigate and exploit MD 
of charged particle pairs and characterize it by a recurrence relation between P„ and P n +i- This assumes connection 
between the collisions of multiplicities n + l with n + l collisions of multiplicity n which can be written in the form 

(n + l)P„+i 

p =9(n). (1) 

The class of MD defined by expression of this type has been considered with a view to stimulated emission and 
cascading in Ref. [l5j . The independent emission of particles represented by Poisson distribution is characterized by 
gin) = c. The constancy of g(n) means that creation of an additional particle is independent of number of other 
particles. The stimulated emission of identical bosons obeying Bose-Einstein statistics follows geometric distribution 
for which one has g{n) = c(n + 1 ). This means that emission probability of a boson is enhanced by a factor n + l 
when n bosons are already present in the system. Both examples are special cases of NBD given by the formula 

_ {n + k-l)\ k 
n n\{k-l)\ Q [ q) ■ ( ' 

The distribution depends on two parameters q =< n > /(< n > +k) and k for which the recurrence relation ([TJ) is 
given by the linear dependence 

g(n) = (k + n)q. (3) 

For fixed < n >, the Poisson distribution is recovered with q, k^ 1 — > 0. The geometric distribution corresponds to 
k = 1. Chew et al. proposed [16| a generalization of NBD which gives good description of MD of charged particles in 
e + e _ annihilation [17]. The generalized multiplicity distribution (GMD) has the form 

_ (n + fc-1)! n -k' (1 }k+fc' m 

n ~ (n - k')\{k - l)\ q 1 Q) • ( ) 

The distribution is function of three parameters q — (< n > —&')/(< n > +k), k, and k! . The expression for g(n) in 
this case reads 

(n + l)(k + n) 

9( n ) = / , , 7TTQ- (5) 
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The formula for GMD reduces to the expression for NBD, when k' — 0. For k — 0, GMD becomes Furry- Yule 
distribution (FYD) proposed by Hwa and Lam [18]. The mentioned distributions are fully determined by the functional 
form of g(n) together with the general normalization condition to unity. The normalization is guaranteed for n > k' — 1 
by the linear dependence of g(n) at large multiplicities with q < 1. In the parton cascade picture of multiple production 
the parameters of these distributions can be related with the rate constants of branching processes in the evolution 
equations for probabilities. 



A. Cascade processes with parton recombination 



In a high energy collision with the creation of a multiparticle state the dynamics of the system can be simulated 
as a parton cascading associated with quark and gluon interactions during the interaction time. The corresponding 
parton cascade processes inspired by the elementary perturbative QCD are of three types: a/ parton-parton collisions, 
b/ branching processes such as quark bremsstrahlung, gluon self interaction or gluon splitting, c/ fusion processes 
e.g. gluon annihilation on a quark. Practically all QCD-based models of many-body production involve some form of 
cascading. The stochastico-physical picture of the multiplicity evolution can be described by Kolmogorov-Chapman 
differential-difference (DD) rate equations [l9[ for the probability P n (t). The continuous parameter t is usually 
interpreted as an ordering parameter, QCD evolution parameter, or time. The rate equations were applied to hadron 
physics many years ago in Ref. [20(. Soon afterwards, the properties of MD were studied by solving DD equations in 
terms of parton cascading by number of people [21] • In particular, all possible QCD vertices were taken into account 
in Ref. [22[. Stationary regime in birth and death processes with recombination (confluence) of two (2 — > 1) and 
three (3 — » 1) gluons has been investigated in Ref. |23j |. 

In this paper we point out and exploit somewhat different strategy. We try to establish minimal number of such 
processes in the cascade picture of particle production which could account for both the general character of MD 
and its complex structure observed at high energies. A successful description of experimental data requires essential 
features of the cascade evolution to be introduced into the phenomenology. We attempt to estimate those features 
which mostly influence the building structure of MD at high energies. This involves processes of parton recombinations 
at the end of the cascade just near breakdown of the confinement. The last phase of the parton evolution is considered 
to have large impact on the form of MD of the produced hadrons. We suppose that in the final stage there exists some 
kind of a stationary regime which occurs at the transition between parton and hadron degrees of freedom. In this 
regime, soft partons intensively exchange their momenta to reach an "momenta uniformity" before their conversion 



into observable particles [23 1 



In order to introduce features of parton recombination into a multiplicity description, we consider a stationary 
regime of DD equations of Kolmogorov-Chapman type. Hereafter we refer to the partons as objects which cascade 
in "time" and regard the terms particle and parton as interchangeable. The binary parton-parton scatterings do not 
change the number of particles during the system evolution but can influence neutralization or screening of the color 
flow. They can act as supporting processes which initiate branching or fusion of the nearby partons. Such connection 
represents multiparton interactions contributing to a change of the particle number by one. The simplest interactions 
of this type are recombinations of two or three particles. Here we consider the process of two parton incremental 
recombination (2 — > 3) together with three parton incremental (3 — > 4) and three parton decremental (3 — > 2) 
recombination processes. The confluence of two and three partons, 2 — > 1 and 3 — > 1, is supposed to be negligible 
with respect to other recombination processes at the end of the cascade. Together with the particle production 
that depends on particles already produced, we consider independent immigration (0 — > 1) and absorption (1 — » 0) as 
another source of partons influenced by bulk properties of the system created in the collision. The accumulated energy 
is transformed into partons in this manner, or on the contrary a parton can be melted and absorbed by the expanding 
system again. In the case of independence of each cascade process, the corresponding DD evolution equations for the 
probability P n {t) read 

P n = + aoPn-i - a P n 

+ /3 (n+l)P n+ i - p nP n 

+ a 2 (n-l)(n-2)P n _i - a 2 n(n-l)P„ 

+ o 3 (n-l)(n-2)(n-3)P n _ 1 - a 3 n(n-l)(n-2)P n 

+ j 8 2 (n+l)n(n-l)P n+ i- j 9 a n(n-l)(n-2)P n . (6) 



The coefficients ao, a 2 , 03 are the corresponding rates for the processes — > 1, 2 — >• 3, 3 — > 4, respectively. Similar 
parameters for degradation of the particle number in the processes 1—^0 and 3 — >■ 2 are denoted as /?o and /3 2 . Using 
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the definition of the generating function 

oo 

Q(w) = J2PnW n , (7) 

we rewrite the system of DD equations ([6]) in terms of Q(w). The corresponding stationary solution (P n = 0) satisfies 
the third order differential equation 

^/ x „ dQ(w) ?d 2 Q(w) , o „.d 3 Q(w) 

a Q(w) - Po-^- 1 + a 2 w 2 J\ ' + (a 3 w 3 - (5 2 w 2 ) J\ ' = 0. (8) 
aw aw z aw 3 

A regular solution of this equation is proportional to the generalized hypergeometric function 

Q(w) = 3 F 2 [ a 1 ,a 2 ,a 3 ;bi,b 2 .—w ) (9) 
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The complex constants at and bj are given in terms of the real parameters uq, a 2 , a^, /3o, and j3 2 (see Appendix A). 
The ratio ([1]) is given as follows 

(n + l)P n +i _ ag+ a 2 n(n - 1) + a 3 n(n - l)(n - 2) 

This recurrence relation together with the normalization condition for P n fully determines the formula for MD which 
we refer to as the generalized hypergeometric distribution (GHD). The distribution depends on four parameters ao/ '/3 2 , 
a 2 /(3 2 , a.3/ /3 2 , and fta/P 2 which are ratios of the rate constants for the corresponding cascade processes. 

In the next section we exploit GHD to describe MD of the charged particles produced in pp/pp and e + e~ interactions. 
The formula (|10[) has various limits applicable for both reactions at different energies and phase space regions. If 
ao -> and j3o — > simultaneously, GHD reduces to NBD. In that case g(n) is given by ([3]) with the parameters 
q = a 3 /f3 2 and k = a 2 /a 3 — 2. On the contrary, if ao and j3o are both large enough relative to other parameters, 
a Poisson-like peak appears in MD at small n for which g(n) ~ ao//3o- There exists a region where GHD can be 
narrower than Poisson distribution. This happens at low energies where the rate constants a 2 and a 3 are small or 
even vanish. 



III. ANALYSIS OF DATA 



The statistics of the charged particle MD is governed by the requirements of conservation laws in each collision. 
Due to the global charge conservation, there are only even multiplicities when dealing with experimental data in the 
full phase-space. On the other hand the charge conservation cannot be satisfied by the stochastic processes which are 
independent of the charge. In particular, GHD, as defined for all non-negative integer values of n, cannot directly 
represent the full phase-space data on multiplicities which are always even. There are two main conceptions how to 
tackle the problem. According to the often used procedure 0] , the probabilities for only even integers n are taken from 
a theoretical distribution, then are renormalized and compared to the measured values of P n . The second approach 
assumes to deal with MD of particle pairs and compare it to a theoretical distribution for all values of n. Both 
methods give different values of the parameters entering the distribution and differ also in capability of description of 
experimental data. Here we use the second method and apply GHD and NBD (for comparison) to the distribution of 
particle pairs. 

We have analysed experimental data on MD of charged particles produced in the non-single-diffractive (NSD) 
pp/pp collisions in the full phase-space. The analysis was performed with the high energy data measured by the E735 
Collaboration [10] at ^/s = 1800, 1000,546, and 300 GeV, by the UA5 Collaboration @at^ = 900 and 200 GeV, 
and by the ABCDHW Collaboration Q at ^ = 63, 53, 45, and 30 GeV. The study includes also the FNAL fixed 
target data [2(1-29] at y/s = 27.6, 23.8, 19.7, and 13.8 GeV corrected for diffraction cross section s fgoll . the Serpukhov 
data at ^fs = 11.5 GeV [3l| with diffraction corrections [32| and the low energy CERN data |33| measured by the 
BHM Collaboration at y/s = 6.8 and 4.9 GeV. The results of the analysis with the CERN-MINUIT program are 
shown in Table I and in Figs. 1-2. As can be seen from Table I, the description of MD of particle pairs by NBD is 
in most cases unsatisfactory including the data from ISR and lower energies. Though MD of the charged particles 
can be approximated by NBD up to the energy y/s = 200 GeV pretty well, NBD is not able to account fully for the 
narrower distribution of the particle pairs. At the energies lower than y/s ~ 20 GeV, the distribution becomes even 
narrower than Poisson distribution. This corresponds to negative values of the parameters k and q, meaning that NBD 
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TABLE I: The results of analysis of MD of the charged particle pairs in the NSD pp/pp collisions in the full phase-space. 
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transforms to a binomial one. By virtue of its four parameters, GHD describes the complex structure of the charged 
particle MD emerging in the high energy pp/pp collisions in the full phase-space sufficiently well. This is illustrated 
in Fig. 1(a) on data [10] from the E735 Collaboration where a peak at low n and a shoulder at large multiplicities is 
visible. The evolution of the observed structure with the energy y/s is shown in more detail in Fig. 1(b). Here the 
relative residues of MD with respect to the NBD parametrization are depicted. The residues are mutually shifted by 
factors 2 for single energies. The lines correspond to the description of data with GHD. As one can see from Fig. 1(b), 
the shoulder broadens with the energy y/s and its maximum moves towards larger multiplicities. This corresponds 
to an increase of the parameter which is ratio of the rate constants for the processes 3 — > 4 and 3 — >■ 2. The 

measurements of MD by the UA5 collaboration at yfs = 200 and 900 GeV in the full phase-space lay systematically 
below the E735 data [10| for large n. The discrepancy has a consequence that the parameter a^/ ' fii is negative for 
the UA5 data. Therefore we set it to null in our analysis in this particular case (see Table I) . 

The description of data on MD by GHD in the NSD proton-proton collisions at the ISR energies is shown in 
Fig. 2(a). The relative residues with respect to the NBD parametrizations are depicted in Fig. 2(b). The residues are 
mutually shifted by factor 2 for single energies. One can see from Fig. 2(b) that NBD does not represent accurate 
parametrization of MD of particle pairs even at the ISR energies. Especially for low n the residues od data relative to 
NBD show remnants of the peaky structure which is clearly visible in the TeV energy region. The solid lines represent 
the parametrization of data by GHD. The small experimental errors allow good determination of 02/^2 and 0:3/ P2 
in this region. The values of the parameters are non-zero at the ISR energies but they are smaller than at the TeV 
energies. Similarly, the parameters a a / ' P2 and P0/P2 are non-zero though both are relatively small. This is once again 
a reformulation of the statement that NBD is not sufficient to describe MD of the particle pairs nevertheless at the 
ISR energies it provides much better approximation of data than in the TeV energy region. 

The analysis of MD below yfs ~ 20 GeV showed that the parameter 013 becomes negative. Therefore we set it equal 
to null in this region. It means that the process 3—7-4 dies out as the first at low energies. For still smaller a/s there 
is not enough energy even for the process 2 — >■ 3 and GHD becomes two parameter distribution. In this region, MD 
of particle pairs is extremely narrow. The non-zero value of 162 makes it narrower than the Piosson distribution. This 
means that the parton recombination process 3—^2 remains active till the very small energies though it brings along 
only the diminution of the particle number. 
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FIG. 1: (a) MD of charged particles in the NSD pp collisions at energies i/s = 300 — 1800 GeV. The data [T3| are multiplied 
by the indicated factors, (b) The relative residues of MD with respect to the NBD parametrization. The residues are mutually 
shifted by the factor of 2 at different energies. The lines represent description of the data by GHD. 




FIG. 2: (a) MD of charged particles in the NSD pp collisions at the ISR energies. The data [25|] are multiplied by the indicated 
factors, (b) The relative residues of MD with respect to the NBD parametrization. The residues are mutually shifted by the 
factor of 2 at different energies. The lines represent description of the data by GHD. 
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A. Structure of multiplicity in limited phase-space regions 



Recently, the CMS Collaboration provided results of systematic measurements of charged particle multiplicities in 
pp collisions at the LHC. The data Q were accumulated in five pseudorapidity ranges from \rj\ < 0.5 to < 2.4 at 
the collision energies y/s = 900, 2360, and 7000 GeV. The measurements of MD in the restricted phase-space regions 
can serve as a sensitive probe of the underlying dynamics in various phenomenological models. 

For an adequate description one needs to make additional assumptions when projecting the full phase-space data 
onto the smaller pseudorapidity windows \t]\ < r] c . The odd-even effect of the charged particle distribution P^ 1 is 
smeared out with the decreasing ij c and the distribution fills all values of n. We avoid the complication concerning 
the charged particle measurements in the limited phase-space regions and collect the data in the neighbouring even 
and odd bins as follows 



P n = P 



ch 



ch 



1,2... 



(11) 



so simulating the distribution of particle pairs. It would be more correct to deal with MD of negative particles which 
is essentially the distribution of the charged particle pairs. Nevertheless, application of GHD to the distribution (fTT|) 
allows us to establish main trends which characterize MD of negative particles. Here we study the dependence of the 
parameters ao//?2i ^2/^2, an d on the srze of the pseudorapidity span < r/ c . 

In Fig. 3(a) we show MD of charged particles in the limited phase-space regions measured by the CMS Collaboration 
in pp collisions at y/s — 7000 GeV. The CMS data for the five pseudorapidity ranges from \rj\ < 0.5 to \rj\ < 2.4 are 
collected in the neighbouring even and odd bins in multiplicity (jllj) and multiplied by the powers of 0.2 for different 
rj c . The data confirm existence of the complex structure observed in MD by the E735 and UA5 Collaborations at 
energies y/s = 900-1800 GeV in the full phase-space. At y/s — 7000 GeV, the peaky form at low n and a shoulder at 
high multiplicities is seen in the all five measured ?7 c -windows. The peak becomes most prominent and the shoulder 
most wide for \rj\ < 2.4. The evolution of the structure with r] c is demonstrated in Fig. 3(b) in more detail where 
the relative residues of the data with respect to the NBD parametrization are depicted. Fore the sake of clarity the 
residues are mutually shifted by unity for different pseudorapidity intervals. The solid lines represent a description 
of the data by GHD. One can see from Fig. 3(b) that at the energy y/s — 7000 GeV the residual structure survives 
even down to rj c = 0.5. As a consequence, NBD is not sufficient to describe the data well enough neither for small 
pseudorapidity windows at this super high energy. 
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FIG. 3: (a) MD of charged particles in the limited phase-space regions from the NSD pp collisions at y/s = 7000 GeV. The 
data H are accumulated in the neighbouring even and odd bins and multiplied by the indicated factors, (b) The corresponding 
relative residues with respect to the NBD parametrization. The residues are shifted mutually by unity for different rj c . The 
lines represent description of the data with GHD. 
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FIG. 4: (a) MD of charged particles in the limited phase-space regions from the NSD pp collisions at \fs = 2360 GeV. The 
data H are accumulated in the neighbouring even and odd bins and multiplied by the indicated factors, (b) The corresponding 
relative residues with respect to the NBD parametrization. The residues are shifted by unity for different r/ c . The lines represent 
description of the data with GHD. 
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FIG. 5: (a) MD of charged particles in the limited phase-space regions from the NSD pp collisions at \/s = 900 GeV. The 
data [3] are accumulated in the neighbouring even and odd bins and multiplied by the indicated factors, (b) The corresponding 
relative residues with respect to the NBD parametrization. The residues are shifted by unity for different r] c . The lines represent 
description of the data with GHD. 



9 



The CMS data measured at the energies \/s = 2360 GeV and y/s = 900 GeV are presented in the same fashion 
in Fig. 4 and Fig. 5, respectively. A similar structure as in Fig. 3 is visible also at these energies. The systematic 
measurements of the CMS Collaboration allow us to study the evolution of the observed structure with y/s. One can 
see from Figs. 3(b), 4(b), and 5(b) that the peak at low n becomes more distinct with the increasing energy. 
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FIG. 6: (a) MD of charged particles obtained by the ALICE Collaboration [II] at the energies y/s = 900, 2360, and 7000 GeV 
in the central pseudorapidity region \r}\ < 1 in the sample n c h > 0. (b) MD of charged particles measured by the ATLAS 
Collaboration [141 ] at the energies yfs = 900 and 7000 GeV in the pseudorapidity region < 2.5 with the selection pr > 
100 MeV and n c h > 2. Data are accumulated in the neighbouring even and odd bins and multiplied by the indicated factors. 
The lines represent description of the data with GHD. 
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FIG. 7: (a) The relative residues of the data from Fig.6(a) with respect to the NBD parametrization. (b) The relative residues 
of the data from Fig. 6(b) with respect to the NBD parametrization. The residues are shifted by unity for different y/s. The 
lines represent description of the data with GHD. 
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At the same time the shoulder widens and its maximum moves towards larger multiplicity. On the other hand one 
can see form Figs. 4(b) and 5(b) that, for the smallest window \rj\ < 0.5, the multiplicity structure vanishes and 
the residues of MD with respect to the NBD parametrization become flat at both these energies. This is in accord 
with the conclusions of Ref. 0) , namely that at y/s = 900 and 2360 GeV the experimentally measured multiplicity 
distributions are well described by NBD for r\ r = 0.5. 

The ALICE Collaboration presented data [13| on MD of charged particles in pp collisions at the energies y/s = 900, 
2360, and 7000 GeV. Figure 6(a) shows the experimentally measured MD collected in the neighbouring even and odd 
bins in multiplicity (|11[) in the central pseudorapidity region I77I < 1. The data are from an event class where at least 
one charged particle in the measured pseudorapidity range is required. The depicted distributions are multiplied by 
the powers of 0.2 for different y/s. The corresponding relative residues with respect to the NBD parametrization are 
presented in Fig. 7(a). The residues are mutually shifted by unity for different energies. The full lines represent 
description of the ALICE data by GHD. One can see from Fig. 7(a) that data on MD for \rj\ < 1 at the energy 
y/s = 900 GeV can be well described by NBD. At y/s = 2360 GeV, NBD is still good description of data with the 
exception at low n where a peaky structure begins to emerge. The peak at low multiplicity is clearly visible especially 
at y/s — 7000 GeV. One can see from Fig. 7(a) that NBD overestimates experimental data for high multiplicities 
(n > 55) at this super high energy, as was already observed in Ref. [13j. 

The ATLAS Collaboration measured the charged-particle MD [ijj in different phase-space regions for various mul- 
tiplicity cuts at three LHC energies. The most inclusive phase-space region covered by the measurements corresponds 
to the conditions \r/\ < 2.5, pr > 100 MeV and n c h > 2. We analyse the ATLAS data obtained under the above r] c 
and pr selections at y/s = 900 and 7000 GeV. Figure 6(b) shows the experimentally measured MD collected in the 
neighbouring even and odd bins in multiplicity (|11[) . The distribution at y/s = 900 GcV is multiplied by the factor 
of 0.2. Figure 7(b) shows the corresponding relative residues with respect to the NBD parametrization. The residues 
at y/s = 900 GeV are shifted down by the factor of unity. The full lines represent approximation of the ATLAS data 
by GHD. The distinct peak at low multiplicities in the ATLAS data represent stringent criterion for description of 
the whole shape of the distribution. Values of the parameters and processes included in GHD are strongly restricted 
by the peak's position, its width and hight as well as by the form of the shoulder at high n. The four parametric 
GHD accounts for smooth transition from the peak region to the shoulder at high multiplicities. Both structures in 
the ATLAS data are measured at such level of the experimental errors which exclude other processes to be present 
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FIG. 8: (a) MD of charged particles in the full phase-space and in the limited phase-space regions from the NSD pp collisions 
at y/s — 900 GeV. The data f2|| are accumulated in the neighbouring even and odd bins and multiplied by the indicated 
factors, (b) The corresponding relative residues with respect to the NBD parametrization. The residues are shifted mutually 
by the factor of two for different r\ c . The lines represent description of the data by GHD. 
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in GHD but those considered. Beside this, we have checked that application of GHD to the MD of particle pairs is 
substantial because when applied to the distribution of all charged particles, the form of data is not reproduced cor- 
rectly. The measurements performed by the CMS, ALICE, and ATLAS Collaborations show that the peaky structure 
of MD is clearly seen in pp collisions in the limited phase-space regions at the LHC energies. The structure becomes 
more distinct with the increasing width of the pseudorapidity window and is most pronounced at the highest energy 
v/5 = 7000 GeV. 

A similar trend in the pseudorapidity dependence of MD is visible in data [24[ measured by the UA5 Collaboration 
in the NSD pp collisions at y/s = 900 GeV. Figure 8(a) shows the UA5 data collected in the neighbouring even and odd 
bins in multiplicity (| 1 1 1) in the full phase-space and in four smaller pseudorapidity regions. The depicted distributions 
are multiplied by the powers of 0.1 for different rj c . Figure 8(b) shows the corresponding relative residues with respect 
to the NBD parametrization. The residues in the smaller pseudorapidity windows are mutually shifted down by the 
factor of two. The full lines represent description of the UA5 data by GHD. Starting from the window \n\ < 1.5 a 
small peak at low multiplicity emerges in the shape of MD. It evolves with the window size and becomes best visible 
in the full phase-space. As can be seen from Fig. 8(b), the residues with respect to NBD are perfectly flat in the 
smallest window \n\ < 0.5 meaning that NBD describes MD in this region well. 



B. Multiplicity distribution in e + e annihilations 

The study of multihadron production in e + e~ annihilations can provide additional information on correlation 
structures which are encoded in MD in an integrated form. The qualitative picture of MD of charged hadrons in 
leptonic processes differs in many properties form pp/pp interactions. The average multiplicities are higher and the 
distributions are narrower in e + e~ annihilation in comparison with the hadron collisions at the same energies. At 
phenomenological level, a comparison of both types of the interactions may reveal which features reflect difference in 
the initial state and which are intrinsic to the parton-hadron transitions. 

Using formula (TlU|) for GHD, we studied the inclusive samples of MD of charged hadrons produced in e + e~ annihi- 
lations in the full phase-space at various collision energies y/s. Similarly as in the previous section, GHD is applied to 
the distribution of particle pairs for all multiplicities n > no- The value of no is the minimal number of the charged 
particle pairs measured in experiment. The results of the analysis are presented in Table II and Figs. 9-10. We 
have analysed data on MD of charged particles obtained by the OPAL Collaboration at the centre-of-mass energies of 
y/s = 172, 183, and 189 GeV [34j- The data with high statistical precision correspond to the energy region sufficiently 
far beyond Z° peak. The three data samples are measured in the multiplicity range which begins with n c h = 8 i.e. 
with no — 4 particle pairs. The analysis includes also data on MD measured by the OPAL Collaboration at lower 
energies, y/s = 161 GeV [35| . y/s = 133 GeV (36| . and at the energy of Z° peak, y/s = 91 GeV |37|- The application of 
GHD to the OPAL data shows that the parameter /?o//32 is compatible with zero and, therefore, it was set to null in 
this high energy region. In such a case, GHD becomes three parameter distribution for which the recurrence relation 
(TTJ) takes the form 

(n + l)P n+1 a _ _ 

= — ; -r + a 2 + a 3 (n - 2), a t = cii p 2 . (12) 

P n n(n - 1) 

In this limit the distribution reflects relatively sharp increase of P n at low n > ng typical for e + e _ annihilations at 
high energies and accounts for a "negative binomial tail" at large multiplicities. 

The relative residues of MD of charged particles measured by the OPAL Collaboration with respect to the descrip- 
tion of the data by GHD with three parameters ([12]) are depicted in Fig. 9(a). The residues are nearly zero reflecting 
acceptable parametrization of the experimental data. Similar holds for the data [38[ obtained by the ALEPH Col- 
laboration at the energy y/s = 91 GeV (see Table II). As shown in Ref. [13], a good description of the OPAL data 
can be obtained also by GMD characterized by the recurrence factor (|5|). An exception represent data on MD mea- 
sured by the DELPHI Collaboration at the energy of Z° peak in the sense that their description requires a non-zero 
value of the parameter /?o//?2- As seen from Table II, this parameter is required als o by data at lower energies. The 



relative residues of MD of charged particles measured by the DELPHI [39], AMY [40(, TASSO [4l|, and HRS [42 1 
Collaborations with respect to the four parametric GHD are depicted in Fig. 9(b). In all cases a good description is 
obtained. 

We analysed also data 43] on MD in e + e~ annihilations measured by the MARK I Collaboration in the low energy 
region. The application of GHD to the multiplicity data for y/s < 10 GeV gives negative values of the parameters 
oW 7?2 and a 2 / P2 with an over-parametrized description. Therefore we set a 3 = a 2 — at these energies. In that case, 
similarly as for pp collisions at very low energies, MD in e + e~ annihilation can be characterized by the recurrence ([1]) 
with 

g- x {n) = Po+fon(n- 1), ft = ft/a . (13) 
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FIG. 9: The relative residues of MD of charged particles produced in e e annihilation in the full phase-space with resp ect 
to description by GHD. (a) The data measured by the OPAL Collaboration at y^s = 189, 183, 172 GeV are from Ref. [34|. at 
yfs = 161 GeV from Ref. [U, at y^s = 133 GeV from Ref. [Uland at ■ s /a = 91 GeV from Ref. '37]. (b) The data measured 
by the DELPHI Collaboration at ^/s = 91 GeV are from Ref. [3Sj. The combined data measured by the AMY Collaboration 
at s/s = 57 GeV are from Ref. The data measured by the TASSO Collaboration at Js = 43.6, 34.8, and 22 GeV are from 

Ref. [4l[. The data measured by the HRS Collaboration at ^/s — 29 GeV are from Ref. [42] ]. 
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FIG. 10: (a) The relative residues of MD of charged particles produced in e + e~ annihilation in the full phase-space with 
respect to description by GHD. The data measured by the TASSO Collaboration at yfs = 14 GeV are from Ref. [41J. The data 
measured by the MARK I Collaboration at ^/s — 7.4, 4.8, 4.0, and 3.0 GeV are from Ref. (b) The energy dependence of 

MD of charged particles produced in e + e~ annihilation in the full phase-space. The depicted data from the OPAL |34|. DELPHI 
[3§|, AMY [40l| . and HRS [11] Collaborations were measured at the energies yfs = 189, 91, 57, and 22 GeV, respectively. The 
lines represent description of the data by GHD. 



Because of the non-zero values of P2 in (|13p . GHD is narrower than Poisson distribution in this region. Figure 10(a) 
shows the relative residues of MD with respect to GHD at low energies. A summarizing illustration of the description 
of MD of charged particles produced in e + e~ annihilations at some energies is presented in Fig. 10(b). The full lines 
represent GHD with parameters quoted in Table II. 
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TABLE II: The results of analysis of MD of the charged particle pairs in e + e~ annihilations in the full phase-space by GHD. 
The minimal number of the pairs is no. If not quoted, the corresponding parameter of GHD was set to null. The errors 
correspond to the quadratic sum of the statistical and systematic uncertainties of data when both are published. 
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IV. DISCUSSION 

The motivation of the present study of multiple production originates from the experimental observation that in 
pp/pp interactions at high energies a structure in the charged particle MD emerges both in the full phase-space and in 
the limited phase-space regions. We concentrate in obtaining a plausible description of the observed structure which 
is distinctly visible in new data from the LHC in the super high energy domain. The phenomenological analysis is 
based on a scenario of multiparticle production in terms of parton cascade processes. The proposed approach aims 
to grasp some qualitative features of parton to hadron transitions which may be important at the end of the parton 
cascading. The observable shape of MD is assumed to be influenced mostly by the soft particles produced in the final 
stages of the cascade development. Besides the ordinary birth (0 — > 1) and death (1 — » 0) process, we consider the 
multiparton incremental (2 — > 3), (3 — > 4) and decremental (3 — » 2) recombination processes which are supposed to 
contribute significantly to the multiplicity build up. Such kind of two and three parton recombination interactions 
in the final stage of the cascade evolution can be justified by the physical requirements of color neutralization and 
reaching of an approximate "momenta uniformity" at hadronization. 

The phenomenological background of data description relies on DD evolution equations which include the terms 
corresponding to the suggested processes. A stationary solution of the equations gives the recurrence relation (|10p 
defining a distribution which we refer to as GHD. The essential ingredients of the analysis are four parameters of GHD 
which are the ratios cvo//32, 0*2/ @2, a3//?2, and /?o//32 constructed from the rates at and Pi of the incremental and 
decremental parton cascade processes, respectively. The energy dependence of the parameters in the full phase-space 
which follows from the performed analysis is shown in Figs. 11 and 12. The empty symbols represent the values 
obtained from data on MD in pp/pp interactions. The full symbols correspond to e + e~ annihilations. The regular 
behaviour of all four parameters is guaranteed by the non-zero values of P2 ■ It means that the recombination process 
(3 — s- 2) is present in both reactions at all energies. The ratios ao//?2 and Po/fo reflect features of data connected 
with the different initial state in the lepton and hadron collisions. They are proportional to the rates ao and /3o of the 
immigration and death processes governed by the interactions of partons with the rest of the system as a whole. We 
will refer to them as parameters of the I. type. As one can see from Fig. 11(a), the parameter ao/ '$2 increases with 
y/~s for both reactions. This means that an increase in the collision energy results in the relative enhancement of the 
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FIG. 11: The energy dependence of the parameters ao//92 and of GHD obtained from the analysis of data on MD in 

the full phase-space. The full and empty symbols correspond to e + e~ and pp/pp interactions, respectively, (a) Ratio of the 
corresponding rates for the process — > 1 of parton immigration and the process 3 — > 2 of parton recombination, (b) Ratio of 
the corresponding rates for the process 1 — > of parton absorption and the process 3 — > 2 of parton recombination. 
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FIG. 12: The energy dependence of the parameters ai/fii and az/fii of GHD obtained from the analysis of data on MD in 
the full phase-space. The full and empty symbols correspond to e + e~ and pp/pp interactions, respectively, (a) Ratio of the 
recombination rates for the processes 2 — > 3 and 3 — > 2. (b) Ratio of the recombination rates for the processes 3 — > 4 and 
3^2. 



parton immigration (0 — > 1) with respect to the three parton decremental recombination (3 — » 2). Figure 11(b) shows 
similar trend in the parameter Po/ f3% for pp/pp collisions. The relative increase of the rate /3q with yfs points to the 
intensive absorption of partons in the bulk of the expanding system formed in hadron collisions at high energies. A 
different picture is foreseen in e + e~ annihilations. Here the parameter P0/P2 reaches a maximum at y/s ~ 40 GeV and 
beyond the energy of Z° peak it drops to zero (see Table II). This suggests that, in contrast to the hadron collisions, 
the parton absorption becomes negligible in the e + e~ annihilations at high energies. 

Another difference in the behaviour of the parameters of the I. type for the lepton and hadron collisions is in their 
absolute values. In the region below yfs ~ 40 GeV, both parameters ao//?2 and /3o//?2 are larger for e + e~ annihilations 
in comparison with pp/pp interactions. There are indications from Figs. 11 and 12 that, at high energies, this tendency 
may be opposite. A reason for such a behaviour may be influence of jets on the multiplicity structure. While in the 
e + e _ annihilations at low energies the immigration rate arj stemming from two (or few) jets prevails the parton 
immigration in the hadron collisions, the multitude of minijets at high energies would result in much larger arj in 
the pp/pp interactions. This in turn would lead to the higher absorption rate /3rj in the hadron collisions giving the 
partons more probability to be melted conversely in the complex system with many minijets again. 
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As seen from Figs. 11 and 12, there exists a region in which the parameters a§j fa and fa/ fa can be relatively small 
within the errors indicated. The region corresponds to the pp interactions in the energy interval yfs ~ 20 — 60 GeV. 
The small values of ao and fa mean that GHD depends mostly on the parameters oii/ fa and a^/ 'fa. In such case the 
data can be relatively well approximated by NBD except for a few low values of n (see Fig. 2(b)). 

The energy dependence of «2 1 fa and 0:3 /fa is shown in Fig. 12. The ratios depend only on the rates of recombination 
processes which are assumed to be active at the stage of parton-hadron conversions. The parameters characterize 
features of data connected with a breakdown of confinement and onset of hadronization. We denote ct2 /fa and 
0:3/ fa as parameters of II. type. Within the errors indicated, both parameters reveal approximately the same energy 
dependence common for e + e~ and pp/pp collisions. Exceptions make the values for UA5 data [2^ which result from 
discrepancies at high multiplicities pointed out in Ref. [l(| ■ As seen from Tables I, II and Fig. 12(a), the parameter 
ail fa has a threshold in the region yfs ~ 7 GeV. Afterwards it becomes larger than unity and continues in rapid 
growth at high energies. This means that, with increasing yfs, the rate of the recombination process (2 — > 3) prevails 
still more and more the rate of the inverse process (3 — > 2). The parameter 0:3/ fa has a threshold in the region 
yfs ~ 20 GeV. It grows with the energy and reaches the value of 0.6 at y/s ~ 1 TeV. 
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FIG. 13: The pseudorapidity dependence of the parameters and fio/fiz of GHD obtained from analysis of MD in pp/pp 
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3 — > 2 of parton recombination, (b) Ratio of the corresponding rates for the process 1 — > of parton absorption and the process 
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FIG. 14: The pseudorapidity dependence of the parameters ct^f fii and 0:3/^2 of GHD obtained from analysis of MD in pp/pp 
interactions at different y/s. (a) Ratio of the recombination rates for the processes 2^-3 and 3 — > 2. (b) Ratio of the 
recombination rates for the processes 3 — > 4 and 3 — ¥ 2. 
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In contrast to the parameters of I. type, the parameters of II. type reflect features of multiplicity production which 
are common to the e + e~ and pp/pp interactions. The analysis suggests that two and three parton recombination 
processes may be part of an intrinsic property of the parton-hadron transitions. Physical justification for such an idea 
may be connected with the processes of color neutralization. 

Experimental data on MD in the central pseudorapidity windows \r}\ < rj c allow us to study the behaviour of the 
parton processes in the limited phase-space regions. The structure of MD observed in the hadron collisions in the full 
phase-space demonstrate itself distinctly in the limited windows in pseudorapidity if the collision energy is sufficiently 
high. This allows more reliable determination of the ratios of the corresponding rates of single processes in dependence 
on the window size. The fine structure of MD is visible even for small r\ c = 0.5 in pp collisions at yfs = 7000 GeV. 
The data at this energy give strongest restriction on the values of the parameters in the small pseudorapidity range. 
The dependence of the ratios a /(3 2 , a 2 /f3 2 , ot 3 /f3 2 , and P0/P2 on r\ c is depicted in Figs. 13 and 14. The symbols 
represent values of the parameters obtained from the analysis of data measured by the CMS, ALICE, ATLAS, and 
UA5 Collaborations at different energies. For clarity reasons, every figure is divided into four panels. Three panels 
show the values obtained from analysis of data measured at the LHC at \fs = 7000, 2360, and 900 GeV, respectively. 
The values corresponding to the UA5 data are shown on the fourth panel. 

One can see from Fig. 13 that ao/f3 2 and /?o//?2 increase with r\ c at all displayed energies. Both parameters reveal 
weak energy dependence in the depicted rj c region. With the increasing window size, the value of «o becomes still 
larger than /3o- This means that the rate of the immigration process (0 — > 1) grows faster with pseudorapidity than 
the rate of the parton absorption (1 — > 0). Except for the energy y/s — 7000 GeV, the parameters a /(3 2 and (3q/(3 2 
can be relatively small for rj c = 0.5 within the errors indicated. In such a case GHD depends effectively on two 
parameters of the II. type and can be approximated by NBD. This is seen in Figs. 4(b), 5(b), and 8(b) where good 
approximation of data by NBD at -^i = 2360 and 900 GeV for r\ c — 0.5 is demonstrated. Such conclusion is in accord 
with the resume made in Ref. Q. As shown in Fig. 13, the errors of ao/(3 2 and Pq/ j3 2 at y/s — 7000 GeV exclude 
small values of these parameters even for rj c = 0.5. This is why the data cannot be well approximated by NBD. The 
statement rephrases the fact that structure of MD stays beyond the NBD description in this region (see lowest part 
ofFig.3(b)). 

The pseudorapidity dependence of the parameters a 2 / (3 2 and a 3 / '/3 2 is shown in Fig. 14. The ratio a 2 //3 2 reveals 
growing tendency with r\ c and flattens at \fs = 7000 GeV. This observation suggests importance of the two parton 
incremental recombination process (2 — > 3) in the small pseudorapidity windows at this ultra high energy. On the 
contrary, three parton incremental recombination (3 — > 4) falls for small ry c , as shown in the upper left panel in 
Fig. 14(b). The growing tendency of a 3 //3 2 with r\ c is seen in the LHC data at all quoted energies. 

V. CONCLUSIONS 

Charged particle multiplicity distributions in pp/pp collisions have been studied including new data from the LHC. 
The analysis comprises MD in the full phase-space as well in the limited windows in pseudorapidity. At high energies 
the distributions show a relatively narrow peak at small multiplicities and a shoulder in the tail. 

A phenomenological description of the observed structure was proposed. Using techniques based on the solution 
of DD evolution equations relevant for the stochastico-physical picture of particle production a simple formula (|10[) 
for the probabilities of the secondary produced multiplicity has been obtained in a stationary regime. The basic 
ingredients of the scenario are the elementary immigration and absorption of partons and the processes of particle 
recombination. We considered two and there particle incremental (2 — > 3), (3 — > 4) and three particle decremental 
(3 — > 2) recombinations. Physical justification for existence of such kind of processes may be connected with the 
requirements of color neutralization at the end of the parton cascade and reaching of an approximate "momenta 
uniformity" of the soft particles at hadronization. The features such as two and three parton recombination allow to 
change particle number, exchange particle momenta and neutralize color repeatedly just before the conversions into 
the observable hadrons. The corresponding solution of the higher order equation for the generating function based on 
the recombination processes exhibit qualitative properties which are absent in the first order. 

This allowed a quantitative description of the complex structure of data on MD in pp/pp collisions both in the full 
phase-space and in the limited pseudorapidity windows. The phenomenological formula (GHD) was applied to the 
description of the charged particle distributions in e + e~ annihilations at different energies \fs. A good agreement with 
data was obtained. The dependence of the four parameters of GHD on the energy and pseudorapidity was discussed. 
The behaviour of some parameters reveals a universal character which is independent of the reaction type while some 
other parameters depend on it. It was shown that the incremental recombination processes play an increasingly large 
role in the multiplicity production as the collision energy increases. 

Within the approach used and on the basis of the studied material we conclude that data on MD indicate existence 
of certain type of recombination processes correlating particle degrees of freedom which manifest itself at high energies. 
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Appendix A 

The system of the evolution rate equations ([6]) can be solved analytically in a stationary regime. The stationary 
solution for the generation function satisfies the differential equation ([5]). Introducing the substitution 



a 3 
Z = T w 

P2 



the equation for Q(w) — F(z) becomes 



/ 3_ 2,d 3 F(z) i a 2 ^ 2 d?F(z) /3 dF{z) ( a 
[z z ) 



dz 3 



03 



az £ p 2 dz a 3 



It can be rewritten into the known form 



z— + ai 

az 



dz 



do \ \ Z 



<,,) - \z- + bl -l 



z— + b 2 -l 

dz 



F(z) = 



(14) 



(15) 



(16) 



for the generalized hypergeometric function 3^2(01, a 2 , a 3 ; b%, b 2 ; z). Exploiting the decomposition of 3 F 2 {z) into the 
power series of its argument the recurrence relation ([1]) can be written as follows 



(n + l)P n +i _ (oq + n)(a 2 + n)(a 3 + n) q 3 
K ~ (61 +n)(b 2 + n) /V 

Here we used P^|) and the standard relation for calculating the probabilities from the generating function, 

1 d n Q 



Pn = 



n! dw n 



(17) 



(18) 



tu=0 



Now, the recurrence relation (|17l) remains to be expresses in terms of the rate parameters a, and ft. For this purpose 
we perform the differentiations in (|16[) and obtain 



d 2 F(z) 

-j^- + z [(3+ ai +a 2 + a 3 )z - (I + 61 + 62)] ~A 



+ [(l + ai + a 2 + a 3 + aia 2 + aia3 + a 2 a 3 )z — bib 2 ] 



dF{z) 
dz 



a\a 2 a 3 F{z) = 0. 



(19) 



The comparison of the coefficients in the equations (fTS")) and (|T9"| gives 



1 + a\ + a 2 + a 3 + a\a 2 + a\a 3 + a 2 a 3 = 0, 

a 2 a o 
3 + ai + a 2 + a 3 = — , ma 2 a 3 = — (20) 
a 3 a 3 

and 

l + 6!+6 2 = 0, b lb2 = ^. (21) 

P2 

The first system of equations represents the correspondence between the sets of parameters Oi and on. The second 
system connects bi and ft. After performing multiplication on the right hand side of (|17j) and exploiting relations 
([2"U)) and (|2"Tj) . one obtains directly the recurrence relation (|TU)) . 

In the analysis of data we deal with the number of particle pairs n. The minimal number of the pairs no is always 
greater than null. For the NSD pp/pp collisions we take no = 1 in all cases. The values of no for e + e~ annihilations are 
shown in Table II. The recurrence relation (jlOl) begins always with n and the obtained values of P n are renormalized. 
The corresponding generating function reads Q{w) — cw n ° 3 F 2 (z). 
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Appendix B 



TABLE III: The results of analysis of MD of the charged particle pairs in pp collisions in the limited phase-space regions |7j| < r\ c 
by the negative binomial distribution (NBD) and generalized hypergeometric distribution (GHD). The data measured by the 
CMS Collaboration are taken from Ref. [3j. 





Vc 


NBD 






GHD 










(GeV) 




k 


q 


x 2 / NDF 


Qo/^2 




a 3 //32 


/V/3 2 


X 2 /NDF 


7000 


2.4 


1.56 ±0.03 


0.910 ±0.001 


89.4/67 


(2.7±0.6)10 2 


6.8 ±0.6 


0.85 ±0.01 


58 ± 11 


6.4/65 




2.0 


1.66 ±0.03 


0.888 ± 0.001 


140.9/61 


(2.3±0.5)10 2 


7.1 ±0.6 


0.82 ±0.01 


56 ± 10 


7.4/59 




1.5 


1.72 ±0.03 


0.851 ±0.002 


143.5/54 


(1.5±0.4)10 2 


7.0 ±0.6 


0.76 ±0.01 


44 ± 9 


3.9/52 




1.0 


1.68 ±0.04 


0.793 ± 0.003 


92.4/38 


76 ± 24 


6.4 ±0.7 


0.67 ±0.02 


29 ±8 


1.9/36 




0.5 


1.69 ±0.07 


0.652 ±0.006 


51/20 


25 ± 12 


5.2 ±0.9 


0.47 ± 0.05 


15 ±6 


1.7/18 


2360 


2.4 


1.99 ±0.08 


0.858 ± 0.004 


30.9/40 


(3.6 ± 1.7)10 2 


9.6 ±2.5 


0.71 ± 0.06 


83 ±37 


9.8/38 




2.0 


1.89 ±0.09 


0.840 ± 0.005 


22.2/33 


(1.7±0.9)10 2 


7.2 ± 1.8 


0.73 ±0.06 


45 ±23 


8.6/31 




1.5 


1.90 ±0.09 


0.797 ±0.005 


17.7/27 


86 ±50 


6.1 ± 1.4 


0.69 ±0.05 


28 ± 15 


4.4/25 




1.0 


2.2 ±0.1 


0.691 ±0.007 


33.4/22 


81 ±43 


7.0 ± 1.6 


0.51 ±0.07 


32 ± 16 


7.2/20 




0.5 


1.9 ±0.2 


0.55 ± 0.02 


5.8/11 


9± 14 


3.2 ± 1.4 


0.44 ±0.14 


5±9 


4.6/9 


900 


2.4 


2.45 ±0.09 


0.794 ±0.005 


28/34 


(2.8 ± 1.3)10 2 


9.4 ± 2.3 


0.63 ±0.06 


70 ±31 


6.4/32 




2.0 


2.27 ±0.09 


0.769 ± 0.006 


18.7/31 


(1.7±0.9)10 2 


7.8 ±2.0 


0.62 ±0.07 


47 ± 24 


2.9/29 




1.5 


2.4 ±0.1 


0.709 ± 0.007 


23.8/26 


(1.5±0.8)10 2 


8.6 ±2.4 


0.48 ±0.09 


51 ± 25 


1.7/24 




1.0 


2.6 ±0.1 


0.603 ± 0.009 


21.6/20 


72 ±46 


6.8 ±2.0 


0.37 ±0.10 


31 ± 18 


2.7/18 




0.5 


2.4 ±0.3 


0.44 ±0.02 


1.7/10 


6± 16 


2.6 ± 1.6 


0.36 ±0.17 


4± 10 


1.1/8 



TABLE IV: The results of analysis of MD of the charged particle pairs in pp collisions in the limited phase-space regions |t;| < r\ c 
by the negative binomial distribution (NBD) and generalized hypergeometric distribution (GHD). The data measured by the 
ALICE Collaboration are taken from Ref. ■ The NSD data sample is used in the analysis at yfs — 2360 and 900 GeV. 





Vc 


NBD 






GHD 










(GeV) 




k 


q 


X 2 /NDF 




«2/^2 


C*3//32 




X 2 /NDF 


7000 


1.0 


1.18 ±0.02 


0.827 ±0.002 


116.4/37 


29 ±6 


4.7 ±0.3 


0.73 ±0.02 


13 ±2 


9.8/35 


2360 


1.3 


1.75 ±0.07 


0.782 ±0.006 


42.7/29 


(1.6 ± 1.1)10 2 


8.3 ±3.4 


0.59 ±0.12 


53 ±35 


23.4/27 




1.0 


1.67 ±0.08 


0.742 ± 0.008 


19.6/23 


44 ±45 


4.9 ±2.2 


0.65 ±0.10 


18 ±18 


5.3/21 




0.5 


1.66 ±0.14 


0.576 ±0.014 


4.9/14 


12 ±20 


3.3 ± 1.9 


0.46 ±0.17 


7± 12 


3.4/12 


900 


1.3 


2.15 ±0.08 


0.702 ± 0.006 


29.9/25 


(2.1 ± 1.6)10 2 


10.2 ±5.1 


0.39 ± 0.20 


70 ±54 


16.5/23 




1.0 


2.00 ±0.09 


0.656 ± 0.008 


12.3/19 


68 ±64 


6.0 ±2.9 


0.47 ±0.15 


28 ±26 


6.8/17 




0.5 


2.14 ±0.17 


0.458 ±0.014 


3.9/11 


29 ±37 


4.8 ±3.3 


0.15 ±0.32 


19 ±24 


0.8/9 



TABLE V: The results of analysis of MD of the charged particle pairs in pp collisions in the limited phase-space regions \r]\ < T] c 
by the negative binomial distribution (NBD) and generalized hypergeometric distribution (GHD). The data measured by the 
ATLAS Collaboration are from Ref. ,14]. The sample with the selection pr > 100 MeV and n c h > 2 is used in the analysis. 





Vc 


NBD 






GHD 










(GeV) 




k 


q 


X 2 /NDF 


a //32 


02//32 


03/^2 




X 2 /NDF 


7000 


2.5 


1.152 ±0.003 


0.9263 ± 0.0003 


4397/54 


(2.4±0.2)10 2 


7.3 ±0.3 


0.826 ± 0.006 


59 ±4 


45.3/52 


900 


2.5 


2.32 ±0.02 


0.790 ± 0.002 


787/35 


(3.7±0.7)10 2 


10.6 ±1.3 


0.60 ± 0.04 


94 ± 17 


29.2/33 
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TABLE VI: The results of analysis of MD of the charged particle pairs in pp collisions in the limited phase-space regions \rj\ < T] c 
by the negative binomial distribution (NBD) and generalized hypergeometric distribution (GHD). The data measured by the 
UA5 Collaboration are taken from Ref. [24| . 



/ — 




NBD 






GHD 










(GeV) 




k 


q 


X /NDF 


Qo/P2 


/ a 
Q-2/P2 


0-3/P2 


P0/P2 


2 /TVTTiTl 

X /NDF 


900 


5.0 


4.26 ±0.13 


0.799 ± 0.006 


95.2/53 


(5.5±4.8)10 3 


39 ±29 


0.3 ±0.4 


(7.7±6.5)10 2 


5.8/51 




3.0 


2.49 ±0.11 


0.819 ±0.006 


25/42 


(7.1±5.9)10 2 


13 ±7 


0.6 ±0.2 


(1.4 ± 1.2)10 2 


5.6/40 




1.5 


2.14 ±0.13 


0.725 ±0.012 


13.2/26 


(1.3 ± 1.5)10 2 


7.5 ±5.1 


0.56 ±0.21 


(4.4±5.0)10 1 


0.6/24 




0.5 


1.78 ± 0.27 


0.48 ± 0.03 


0.7/10 


(0.2 ± 1.3)10" 4 


1.77 ±0.09 


0.50 ±0.04 


(1.3±9.3)10 -5 


0.3/8 


200 


5.0 


5.52 ±0.37 


0.654 ±0.015 


7.2/29 


(2.9±0.6)10 3 


31.5 ±2.5 




(4.9±0.9)10 2 


0.8/28 




3.0 


4.06 ±0.27 


0.663 ±0.015 


8.7/26 


(1.6±0.4)10 3 


27.3 ±2.5 




(3.5±0.7)10 2 


0.9/25 




1.5 


3.46 ±0.32 


0.54 ±0.02 


3.4/15 


(2.8±0.9)10 2 


13.0 ± 1.6 




(1.0±0.3)10 2 


1.7/14 




0.5 


5.5 ±2.1 


0.19 ±0.05 


0.4/5 


13 ±32 


2.4 ± 1.3 




11 ± 25 


0.7/4 
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